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Abstract 

A well known method for convergence acceleration of a continued fraction K(a n /6„) 
is based on the use of modified approximants S n (ui n ) in place of the classical ones S n (0), 
where u> n are close to tails of the continued fraction. Recently (Numer. Algorithms 41 
(2006), 297-317), the author proposed an iterative method producing tail approximations 
whose asymptotic expansion's accuracy is improved in each step. This method can be 
applied to continued fractions K(a„/6„), where, for sufficiently large n, a n and b n are 
polynomials in n (dega„ = 2, deg&„ < 1). The purpose of this paper is to extend this idea 
to the class of continued fractions K(a n /b n + a' n /b' n ), where a n , a' n , b n , b' n are polynomials 
in n (dega„ = dega^,deg&„ = degb' n ). We give examples involving continued fraction 
expansions of some mathematical constants, as well as elementary and special functions. 



1 Introduction 



Many valued mathematical constants and elementary or special functions have continued frac- 
tion expansions. For example, the following continued fraction (CF, in short) 
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(1.1) 



is very useful expansion of the well-known inverse trigonometric function. If, however, x ~ 1, 
then it converges very slowly and application of some methods of convergence acceleration is 
required (see p. 19, Eq. (11)]). If x ^ 1, then CF (1.1) can be transformed equivalently 
into limit 1-periodic CF of loxodromic type, and one can use the improvement machine — the 
method established by Lorentzen and Waadeland (see [2], [U p. 227]). But if x = 1, then ( |1.1[ ) 
is of elliptic type and, according to [2], one can use its even part and transform it to the CF of 



parabolic type, of the form K(c n /1), where c n + 
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one can use modified approximants S n (— 1/2) — the fixed point method 



+ O (n- 5 ) . Since c n + \ = O (n" 4 ) , 



the convergence of K(c n /1) (see [U p. 219]). However, if c r , 
usually does not succeed. 



+ 



On 



-2 



in order to accelerate 
i , then such a method 



One can verify that CF (1.1) can be effectively computed using the 
method introduced in [3], where the convergence acceleration of CF G C20 U C21 was proposed; 
mCki means the class of CFs of the form K(a n /6 n ) with a n and b n being polynomials in n 
(dego n = k, degb n = I). In this paper, we extend this result to the class of so-called 
two-variant CFs (see 0) of the form 
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(1.2) 
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where a n , a' n and b n , b' n are polynomials in n, such that deg a n = deg a' n = k, deg b n = deg b' n = I. 
One can find a lot of examples of such CFs in the Perron's book [7\. More CF expansions can 
be found in (chronological order): [9], [3] and pp. The authors of pQ have developed the website 
which shows the applications of continued fractions to approximation of many special functions 
(see http : //www . cf hbli ve . ua . ac . be/| ) • 

For example, the following CF expansion holds 



e _fcE cn(t; k) dt 
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(0 < k < 1, x > 0), (1.3) 



where cn(i; k) is the Jacobi elliptic function, very often applied in electro- and magnetostatics, 
fluid dynamics, as well as nonlinear integrable equations (see p. 220, Eq. (15)] or [H p. 283]). 
Although the equation ( |1.3[ ) is known for almost 100 years (the first edition of Perron's book 
[7] was issued in 1913), current computer algebra systems do not use it. For instance, the 
Maple system is not able to compute the integral in equation (1.3) with x = 0.8 and k = 0.9; 
if we replace the upper limit oo by 500, then the correct result is obtained after 1.5 minutes 
of waiting. Thus, the problems which require computing a lot of such integrals are almost 
insolvable using the current mathematical software like Maple or Mathematica. The motivation 
for using aforementioned techniques of convergence acceleration is, among others, the fact that 
CF expansion (1.3) allows us to compute such integrals in a split of a second. 

Our main result is a new method for convergence acceleration of two-variant CFs. Below, 
we give some notation used in this paper. In Section [2j we describe the class T> of considered 
CFs (see ( |2.4[ )). We analyse tails of CFs and try to find out their asymptotic form. Our 
investigations require to study many cases of CFs' subclasses. For the sake of clarity, we moved 
some technical results to |Appendix A We show that derivation of the coefficients of asymptotic 
expansion of tails is quite costly. Thus, we propose an iterative method that improves the 
accuracy of tails' approximations. In Section [3j we describe our method in detail. The main 
result is in Theorem 3.1. Next, in Section [4] we use this result repeatedly in order to obtain 
the algorithm, which is of iterative character, producing tail approximations whose asymptotic 
expansion's accuracy is improved in each step. This idea is summarized in the Algorithm 4.1 



Finally, in Section [5j we report the results of experiments involving some CFs belonging to the 
class V. 



Some notation 



Let S n (oj) denote the n-th modified approximant of CF (1.2); see [3j pp. 7, 56] or [U (1.2.1), 
p. 5] for more details. Clearly, we have 
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We use the symbol to denote the n-th tail, i.e., 
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2 Asymptotics of tails 



Our method for convergence acceleration is based on the obvious fact: if CF is convergent to 
V, then all its tails are also convergent, and 



Sn(f ( 



V. 



(2.1) 



Since S n (uj) is the linear fractional transformation in u, it is of great importance to find the 
best approximation oj n of the tail f^. Then the modified approximant S n (uj n ) may be a very 
good approximation of V. 

In Subsection |2.1[ we describe the class T> of considered CFs. Taking into account their 



'two-variant' character (cf. (1.2)), we deal with odd tails only, 



One can verify that tails u n of CF (1.2) satisfy the following recurrence relation: 
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In Subsection 2.2, we analyse the behaviour of tails u n and state their asymptotic expansion 
(see ( |2.7| )). Also, we give formulae for its leading coefficients (see Theorem 2.1). Let us remark 
that derivation of them requires considering all the subclasses in T>, separately. Thus, we move 
part of calculations to |Appendix A 



2.1 Class V 



We consider some subclasses of CFs of the form (1.2), where 



a n+ i=p- 2 n +p-in + p + ... 



- 2 n + q-in + q + ..., 
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(2.3) 



We say that CF (1.2 ) is of class V^i if Pk,p' k , qi, q[ are leading coefficients in the above expansions. 
In this paper, we consider only k € {1,2} and I E {0,1}, since any CF of class T>k+2,i+i is 
equivalent to some of the class T>u- 

Let us fix some notation related to CFs' subclasses considered here. In this paper, we 
consider the CFs belonging to 



V := T>Y U T>f n UX»n U U vf n U V~ 



10 



'20 



'20 



'2i ; 



(2.4) 



see Table [TJ where the characterization of all the subclasses is given. 
2.2 Asymptotic expansion of tails 

It follows from ( |2.2[ ) that tails u n satisfy the equation (see also [H (1.6)]): 

{b' n b n+ i + a n+ \) X n + b' n X n X n+ i — a n X n+ i — a! n b n+ i = 0. 



(2.5) 



Now, we analyse the asymptotic behaviour of its solutions. Observe that recurrence rela- 
tion (2.5) is bilinear and all its coefficients have the asymptotic expansion in powers of n _1 
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Table 1: Considered subclasses of CFs 



Subclass symbol 


Characterization 


U w C U\q 


n •. - r>' and ^0 ^0 a ( n i 
p_i — Pi ana tt 1 00, uj 

P-i 


/ 

T>ta c p io 


\p'-i\ + \v-\\ 


^20 C ^20 


P-2 =P-2 and 2 ^ ( oo,0] , 

P-2 

where a = (3 = p_i - p / _ 1 - p_ 2 , 7 = ~P-2 <7o 




IP-2I / IjP— 2I 


V 2 1 C D21 


The roots xq, x\ of equation x = v ~% 2 satisfy -F^ -10/ — xi 

q' . + : — 9-1 9-1 



(cf. (2.3)). Therefore, its solutions have the asymptotic expansion in powers in 1 / 2 . Moreover, 
this expansion is at most 0(n 2 ): 



X n = r_ 4 n 2 + r_ 3 n 3/2 + r_ 2 n + r_i ra 1/2 + r + n n 1/2 + C(n x ). 



(2.6) 



More precisely, equation (2.5) has two solutions X' n ,X'^ of the form (2.6). According to 
Paszkowski [BJ, the sequence ix n is a solution X n G {A"^, A""}, such that {A n a n+ i/(6 n+ i+A n+ i)} 



is asymptotically smaller in modulus. It means that the problem in solving the equation (2.5) 
reduces to choosing suitable coefficients tj, such that 



u, 



fi/2 



+ n-u ^ /2 - 1/2 + T 2 ^ + r 3 _ M n^' 2 + . . . . 



(2.7) 



Of course, knowing the coefficients Tj, we can give good approximations of the tails u n , and 



good approximations of CF value, as well (cf. (2.1)). In this Subsection, we give formulae 
for some portion of the coefficients T-a,Ti-u, ■ . ., in the case of any subclass in D, named in 
Subsection 12.11 



First, we put the value n < 4 (cf. ( |2.6[ ) and (2.7)). From (2.3), it follows easily that /x < 2 if 
CF ^ T>2o- In the sequel, we allow to be r_ M = 0. We will show that this may happen only in 
the case of T>f Q or T>\\. Therefore, we call r_ M the beginning coefficient (not the leading one) of 



expansion (2.7). According to Lemma A.l, we choose 



1, if CF e Vf , 
4, ifCFGP 2 = , 
2 otherwise, 



and assume that t_ v = if v > fx. The proof of Lemma A.l is rather technical and requires 
considering all the subclasses in T> separately. 

Second, let us look at the beginning coefficient r_^. Our main results are as follows. Lemma 

0, while Theorem 



A.2 



shows that t-u satisfies quadratic equation cut 2 + /3t + 7 



2.1 



specifies 



which root it is, where a, (3 and 7 are given in Table [2] Also, we give there simple formu- 
lae for r_^; in Lemmata |A.ll |A.3| one can find the formulae for the next coefficients of the 



expansion (2.7) 
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Table 2: Coefficients a,/3,j 



Subclass 


a 




7 


VTo 







-90 P-i 




9o 


P-l -P-i 





V n 


q'-i 


9-1 q'-i 





^20 


q'o 


P-l -p'-i -p-2 


-90 P-2 


u 20 


q'o 


P-2 - P'_ 2 





T>21 


q'-i 


P-2 - p'- 2 + 9-19-1 


-9-1 P-2 



Theorem 2.1. T/ie beginning coefficient is given explicitly in the following way. 
1. In the case ofV^Q = !)•' 



/ qop-i 

, cvj^V «o 1 /90P-1 

r-i := sgn K — v / — ; 

P-l / V 9 



(2.8) 



2. In the case ofT>f Q ([/, = 2): 



T-2 ■= < 



P-l - P-l 

9 



if \p'-i\ < |P-i| 
*/ |P-i| > |P-i| 



5. In the case ofT>u (fi = 2): 
4- In the case ofV^ = 2): 



t-2 := 0; 



(2.9) 



T-2 



2a 



(2.10) 



where a,/3 and 7 are given in Table^ 
5. In the case ofT>2Q (yu = 4): 

'0 



T-4 := < 



P-2 ~ P-2 
9 



if \p'-z\ < IP-2I 
|P- 2 | > IP-2I 



6. In the case of V21 (fJ- = 2): we choose r_2 as the first of the roots {xq,xi} of equation 

(2.11) 



(A.8) if 



P'-2 


< 


P'-2 


-, X 


-, a?i 


9-i 




9-i 



Proof. According to Paszkowski [6] , the sequence of odd tails fi 2 ^ 1 ) of CF Oh is such a solu- 
tion of recurrence equation (2.5), for which the product j( 2 «-i)j( 2n ) i s asymptotically smaller 
in modulus. 
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Before we start, let us observe that setting 



T f(2n-l) f (2n) _ a n+l "n 

On+1 + Un+1 



yields, after some simple calculations that use only Lemmata A.l and A. 2 



1. In the case of £> 10 : 



2. In the case of T>f : 



3. In the case of T> 



ii- 



4. In the case of T>. 



20- 



I n = p-i n V n + 

T-l 



In 



V-i " + C(l) if r_ 2 = 0, 
p_i n + 0(1) otherwise; 



In = < 



V-iP-i 
9-1 9-1 



+ 0(n~ 1 / 2 ) ifr_ 



gL 1 g_i n 2 + O(n) otherwise; 



I n =p- 2 n + p_i -p_2 



P-290 
T-2 



n + C(l); 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



5. In the case of X>|o : 



\p'_ 2 n 2 + 0{n) ifr_ 4 = 0, 
I p-2 n 2 + 0(n) otherwise; 



(2.16) 



6. In the case of T>2\ 



I, 



n 2 + 0(n^ 2 ). 



(2.17) 



T-2 + q-i 

Let us observe that the last formula is well-defined. Otherwise, assuming r_2 = —q~\ and using 
(A. 9), we derive p-2q~\ = 0, which is a contradiction. The formulae (2.13), (2.14) and (2.16) 
easily imply the parts 2, 3 and 5 of the thesis. To complete the proof of the parts 1 and 4, 
it is enough to make use of the two obvious facts. First, if z n = an K + bn x + 0(n u ) for any 
K > A > v, then 

\ Zn \ 2 = \ a \ n 2K + 2{^am + %a%b)n K+x + 0(n max{2X ' K+u} ). 

Second, sgnK^ = sgn(Ka!R& + QaQb) holds for any a, b G C (b ^ 0). In order to complete the 
proof of the part 6, let us observe that inequality (2.11) is equivalent to 



\xoQ-l -P-2\ < \ x l Q-i ~P-2\ ■ 
By multiplying both sides by \p-2 ff-i/si-il) we obtain 

1-1 P-2 



g-lP-2 
; P-2 + X P-2 <?-l 



< 



Q-i 



p-2 +x 1 p- 2 q-i 
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Now, let us note that xqX\ 



g-ip 



(equal to the quotient j/a in the case V21; see Tablej^). 



Hence, we have \xq x\ p-2 + xq Q-i P-2 1 < l^o ^iP-2 + x± q-\ p_2 1 , and thus 



XqP~2 


< 


XlP-2 


xq + q-i 




x\ + q-i 



Therefore, the choice t_2 := xq implies that I n has asymptotically smaller absolute value 
(cf. pTf])). □ 



3 Method 

In this Section, we present the method for convergence acceleration. We consider only CFs of 



subclasses in V, given in (2.4). However, in general, our idea can be applied to any CF whose 



sequence of odd (or even) tails u n has asymptotic expansion of the form (2.7). Then the only 
problem is to derive the initial approximations of the tails u n . 

The main objective is to improve, in each step, the accuracy of approximations of tails 
u n . Of course, one can use the technique of Subsection |2.2| to obtain a suitable number of 



coefficients Tj in the expansion ( |2.7| ). This approach was discovered more than half-century 
ago by Wynn in [ID] : see also [H p. 232]. Let us note that in such a case, each coefficient Tj 
can be expressed in terms of all the previous coefficients Tj—i, Tj'-2> • • • j T -^i an d hence the 
computation of them explicitly is quite costly. For example, in the case of CFG C20 U C21, the 
recurrence formula for the sequence of tails of CF was given in [5j Thm. 2.6]. 



Our goal is to obtain better and better approximations of the tails (2.7) without computing 
the coefficients Tj explicitly. We will construct approximations u[ m ' of tails u n , such that 
are of increasing accuracy, for consecutive values of m, in the sense that: 

{u^ +l) - u„)/(ni m) - un) -> when n -> 00. 

, n = 0(n- m / 2 ), then ui m+1) -u n = 0(n~ m / 2 - 6 ), where 6 > is given 



More precisely, if u{, m ' 



in (|3.11|). Our method is of the iterative character, since we express u[ m+1 ' in terms of 



A.l 



and A. 3 it follows that tails 



and u^+i- Similar concept was given by Lorentzen and Waadeland, and called 'improvement 
machine' (see [U p. 227] or [21 §6], for more details). However, their method cannot be applied 
to all the subclasses in T> and, as numerical examples show, is less efficient than ours. 

The following technique is very similar to the one presented in [5], but is more subtle, 
since we deal with more subclasses of CFs and consider only t heir odd tails , whose asymptotic 
expansion can have powers of n -1 / 2 . However, from Lemmata 
u n mostly have expansion in Only in the case of T>\ 

powers of n _1//2 . Numerical examples show that it is hard to accelerate such a class of CFs. 
Step of iteration 

Now, we describe a single step of the method in detail. Let us assume that we have some 
computable approximations u' n of tails u n , in the sense that 



the expansion (2.7) may contain 



S' n := v! n - Un = cn-™' 2 + 0{n^ m+ ^l 2 ), (3.1) 
where m > is an integer number and c 7^ is an unknown constant. Here, the name computable 



means that we can compute them for consecutive values of n G N. The formula (3.1 ) says that 
asymptotic expansions on u n and u' n are identical up to the term n~^ m ~ 1 ^ 2 , and so we say that 
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u' n has the order m/2. Since the approximations u' n are different from u n , they do not satisfy 
the equation (2.5). Thus, let us define new approximations, 



b' + 



(3.2) 



b n +l + <+l 



(cf. ( |2.2[ )), and observe that they are computable, too. They do not give better results than the 
approximations u' n , because one can check that modified approximants S2 n -i(Un)i 5 , 2n+i(^n+i) 
are equal. Nevertheless, we can do a trick, using both of them in order to obtain new approxi- 
mations u n of higher order. The main results are summarized in Theorem 3.1, and the details 
are in its proof. The trick is the same as in [5] . Roughly speaking, we want to find a relationship 
between 5' n and 

5t'= u t- u n- (3-3) 

It can be easily verified that may be expressed by d~' n+1 as follows 



where 



ipr, 



+i' 



a n a n+l 



(3.4) 



(3.5) 



(a n+ i + b' n b n+ i + b' n u' n+1 ) (a n+ i + b' n b n+1 + b' n u n+ x) 
However, the formula (|3.4|) is useless, since ip n is expressed in terms of the tails u n +i, which 



are unknown and rather incomputable. The major concept is to transform S' n+1 into S' n in 
the relation (3.4) and put u' n+1 instead of in (3.5). This is done in order to obtain the 
relationship of the form 

v'X-i>'Jn = 0(n- m /i-<), (3.6) 

for some computable functions <p' n , ip' n and £ > 0. Then the new approximations u n can be 
found a solution of (3.6), provided that the right-hand side of (3.6) is replaced by zero and u n 
(stemmed from 5' n and <5+) by u'^. Clearly, we express u'^ in terms of u' n , n+, ip' n , ijj' n and obtain 



(^-V4)K-^n) = 0(n- m/2 - C ). 



(3.7) 



Theorem 3.1. Let the approximations u' n of tails u n satisfy (3.1) with m such that 

'o, ifCFeV^ uvf , 

m > < 1, if CF£Vf , 
2, otherwise. 



(3. 



Let k+ be given by (3.2), and let auxiliary functions ip' n , (p' n be given by 



a n a n+l 



(a n+ i + b' n b n+l + b' n u' n+l ) 2 ' 
Then the approximations v!' n defined by 

_ <Pn U t ~ j>'n Un 
Pn ~ fn 

u n = 0(n- m / 2 - e ), 



Pn 



1 + 



m 
2n' 



// ._ 



satisfy 



(3.9) 
(3.10) 
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where 



1, ifCF€VY UV^ , 

2, ifCFevf uvf uv 21 , 

4, ifCFeVu- 



(3.11) 



Proof. Roughly speaking, the idea of the proof is to modify the relation (3.4) in order to obtain 
relation (3.6). 



First, we use the following formula relating 5' n+1 and d' n : 
Cl = (! " In) 5 « + 0(n- m / 2 -«), where £ := 



|, ifCFGPfo, 
2, otherwise. 



(3.12) 



This can be verified by substituting n+1 in place of n in (3.1 ) and using simple algebra. Observe 
that number £ is common for all classes in T> except the class T>^ ; it follows from Lemmata A.l| 
and | A. 3 In the case of V= the formula (3.12) is exactly the same as the one given in [5, 
Eq. (3.7)]. 

Second, let us observe some properties of function ip' n . One can check that 



= cn p + 0(n p 1 ^ 2 ) for some c ^ and p := 



2, in the case of T>n, 
0, in all other cases. 



(3.13) 



This is not obvious only in the case of T>2i- Clearly, the coefficient of n 4 in the denominator of 
ip' n is then equal to p_2 + q'_i (7-1 + 4-\ r -2 and it does not vanish, because of the choice of r_2, 
given in Theorem 2.1 Next, using the definition (3.5) yields 



a 'n bn^n+l 



and thus 



where 



(a«+l + K K+l + K «n+l) (On+l + + K <+l) 



^ n - ^ = cn -fV + 0(n~T V-i/2) (c ^ o), 



2 °n+lJ 



1, in the cases of £> 10 , and P2I) 
p' := / 3, in the case of Pn, 

,2, in the cases of V^q and £>|^. 

Now, using the fact that ^ + ^ + p<m + p' (cf. ( |3.8|)) a nd multiplying equation (3.4) by 
ip' n (some simple algebra is needed) yields the relation ( |3.6| ) for ( := £ + p > 0. Further, let 
us observe that approximations it", given in (3.9), can be found as a solution of the following 
equation, derived from (3.6) by replacing right-hand side by zero and u n (stemmed by S' n and 
6+) by <: 



Observe that equation (3.7) follows immediately by subtracting (3.6) from the above formula. 
To complete the proof, let us state the following result: 

<p' n — tp' n = cn~ v + 0\n~ r, ~ 1 l 2 ) , where c 7^ and 77 := 



1 in the case of "D in , 



1, in the case of T>- 



20) 



(3.14) 



0, in all other cases. 
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Hence, the formulae (3.10) and (3.11) hold. Observe that it is important that c 7^ in (3.14), 
since we need to divide the equation (3.7) by cp' n — ij)' n . 

To verify the formula (3.14), let us determine the beginning coefficients, say o~j, of the 
asymptotic expansion of ip' n — ip' n : 

- i>'n = a + crirT 1 / 2 + cr 2 n _1 + 



By Lemmata A.l and A. 3 we know that the coefficients 02.7+1 (j > 0) may be nonzero only in 
the case of T>f . The proof is based on the following facts: 



in the case of T> 10 : 00 = 0, 01 = 2<? ° ^ 1 7^ 0; 
in the case of T>f : 00 = 1 — ^7 7^ 0; 



p-i 



• in the case of T>\\: 00 = 1; 



in the case of P 2 o : 00 = 01 = 0, 02 



2r- 2 go+P-i-P-i+fP- 2 



P-2 



+ (cf. (LOT)); 



• in the case of T>2Q '■ 



-0 = 1-^^0; 



in the case of 2?2i : 00 = 1 



P-2 P-2 



^0. 



Only the last one needs some explanation. Namely, let us put 



z := (p-2 + q-i q'-i + q'-i T- 2 ) 2 - P-2P-2- 



One can verify that 



(2z - A) 2 = A( P -2 + p-2 + q-i q'-i 



1 \2 



where A is discriminant of the quadratic equation given in Lemma |A.2 In the case of subclass 
X>2i, we have obviously A ^ (see Table [l). Suppose that Co = 0. Using the following simple 
algebra: 



A P - 2 p'-2 = {p-2 +P-2 + q-i q'-i? -A = {A(p-2+ P '-2 + q-i q'-if - A 2 ) /A, 
yields z = 0, and thus (2z — A) 2 = A 2 . This implies that P-2P-2 = 0, which is a contradiction. 



□ 



4 Algorithm 

Now, we summarize all the previous results. We give here the algorithm for convergence accel- 
eration of CFs of the class V (see ( |2.4| )). It is almost fully rational in the sense that it makes 
only arithmetical operations. Only in the case of CFs of subclass T>f the initial calculations 
need to compute square roots. In the main loop of the algorithm, only arithmetical operations 
are used, for all the subclasses in V. 

The algorithm is divided into two parts. First, we calculate initial approximations u u q of 
tails u n of CF, such that u n o — u n = (D(n~ m / 2 ), where 

0, in the cases of Rfo an d ^20 > 



m 



1, in the case of T> 1Q , 

2, in other cases 
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(cf. (3.8)). There we use the results given in Section 2.2 The second and main part of algorithm 
is the iteration which improves, in each step, the accuracy of approximations u n j of the tails u n . 
Clearly, using the Theorem 2.1 repeatedly, we obtain the following relation for approximants 



ii 



u ni -u n = 0(n- m / 2 -i e ). 



Algorithm 4.1. Consider the continued fraction ( |1.2[ ) belonging to V (see ( |2.4[ )). 
Initial step. Let us put the initial approximants u n o as follows: 
a) if CFG Vf , then 



UnQ ■= T_l n 1/2 + T 



where 



/ QQP-1 



r_i := sgn 5J- 



9o 



P-l 



qop-i 
So ' 



to := 



2p' Q - 2q' q + p-i - 2p 



b) if CFG Vf , then 
where 



u n0 := r_ 2 n + r , 



r -2 := < 



0, if \p'_ x \ < 

P-l - P-l 



otherwise, 



to ■= < 



P-i go 
P-i - P-i 

P-lQo . P-1+Pq-Po . (P-i- P-i) Qi 



I P-l - p-l 



+ 



+ 



if \rf_ x \ < \p-i\ 
otherwise; 



c) if CFe £>n, then u n o '■= 

d) if CFG V^ , then 

U n 



P-l 



7' ' 
1-1 



-fi + sgn (fc^j^) a/^4o^ 



2a 



n, 



where a = g , (3 = p_i - -p_ 2 , 7 = -<?oP-2; 
e) if CFG T>f , then 



where 



T_4 = < 



it n0 := t_4 ra + r_2 n, 



if |P- 2 | < IP-2I, 



P-2 -P-2 



otherwise, 



0, 



T_2 = < 2p_ 2 + p'_! - p-l <?i(p-2 - P. 



+ 



Q'o 



if |P- 2 | < \P- 
otherwise; 
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f) if CFGP21, then 
where 



U n o ■= T_ 2 n + T , 



TO 



P-2Q0 + P-iQ-i - {q'-xqo + q' q-i - v'-\ + P-i ~ P-2) T -2 - (<?-i + q'o) T -2 



2 q'_ X T- 2 + p-2 - p'-2 + q'-iq-i 

and T-2 is one of two roots {xo,x±} of equation 

q'-i + (p-2 - p'- 2 + q-i q'-i) x - q-\p'-2 = Q , 

which yields less absolute value 



P'-2 

x k 

q-i 



k = 0,1- 



Main loop. Compute the approximations u n j using the following formula: 



u n ,j+i '■= 



<J) _ l (3) 



n > 1, j > 0, 



where 



TIJ 



gn+1 



K + 1. 



1 



m 



+ j6)n 



-1 



Cj) 



(a n+ i + 6^ 6 n+ i + b' n u n+ltj ) 2 
5 Numerical results 

We give examples of some remarkable mathematical constants and special functions having con- 



tinued fraction expansions. We apply our method (see Algorithm 4.1) in order to approximate 
their value. 



Consider a two- variant CF of the form (1.2) belonging to the class T>. Of course, one may 



always approximate its value using the following classic methods: 

1. computing classical approximants S n (0) for sufficiently large n, 

2. applying the methods given in [2] or [5] to even or odd part of CF, 

3. using equivalent transformation in order to obtain limit 2-periodic CF, and using the 
techniques given in [2], 

4. using the recent methods of Paszkowski [6J. 

See [U §5] for more details on numerical computation of CFs. The first idea should be rather 
the last resort in order to compute any CF. Convergence of CF from the class T> may be very 



slow (see Example 5.1). 

The second method requires to derive even or odd part of CF (see, e.g., [H Thm. 2.19-2.20, 
pp. 86-87]). After that, in order to apply the techniques summarized by Lorentzen in [2], 
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some equivalent transformation (see, e.g., [H Thm. 2.14, p. 77]) into limit periodic CF may be 
needed. However, for some subclasses in V, this process leads to the case that the methods 
given in [2] or [5] are not applicable. For example, the 'improvement machine', mentioned in 
Section [3| can be applied only to CFs of the so-called loxodromic type; see, e.g., [21 p. 73] for 
definition of CFs of loxodromic, elliptic, parabolic and identity type. Clearly, one can verify 
that odd (even) part of CF of the subclasses P^jPJq is equivalent to a limit periodic CF of 
parabolic type. Then the question of its convergence acceleration is much more subtle unless 
it is divergent. In the case of limit periodic continued fraction K(c n /d n ) of parabolic type, one 
can always try to use so-called fixed point method (see [2], [U p. 219]), i.e. computing modified 
approximants S n (— 1/2), provided that, without loss of generality, c n — > — 1/4 and d n — > 1. But, 
in order to accelerate its convergence, there are some additional conditions required (see, e.g., 
[21 Eq. (4.10)]). In examples 5.2 and 5.3 we show a class of CFs (belonging to V^ and T>Yq) 
having odd (even) part of parabolic type, for which 'fixed point' method does not work. 

The third method is quite similar to the previous one, since it depends on the type of limit 
periodic CF. Still the 'improvement machine' can be applied only in the case of CF of loxodromic 
type. Other cases need more careful handling. Anyway, all mentioned methods of computing 
limit periodic CFs do not cover the considered class T> of CFs. 

The fourth method, established by Paszkowski [6], is based on analytical transformation of 
tails of CF. It can also be applied only to some subclasses in T>. For example, it does not handle 
the case of CF with tails having the asymptotic expansion in powers of n -1 / 2 , and so it cannot 
be applied to the subclass T>^ . Moreover, it is quite difficult to compare such analytical method 



with Algorithm 4.1, which is rather of numerical character. 

In the sequel, we give a series of examples showing the application of our method (Algo- 
rithm 4.1) in order to accelerate the convergence of CFs belonging to the class T>. We give 



additional remarks on computing their value by using aforementioned methods, and, if it is 
possible, we show a comparison. 

Before we begin, let us fix some notation to measure the accuracy of the method. Consider 
the CF K.(a n /b n ; a' n /b' n ) of the class V, convergent to the value V. Let acc(x) denote number 
of exact significant decimal digits of x in V: 



acc(x) := — log 



10 



X 

V 



We use the following notation: 



acc(5 ri (0)), 5 nj := acc (5 2n -i(^nj)) , 



where u n j are the values computed by Algorithm 4.1 Let us remark that u n j is an approximation 
of the tail u n = /( 2n_1 ) of CF, and that is why we analyse the modified approximant S2n-i(u n j) 
(cf. (2.1)). We call the number &cc(S2 n -i(u n j)) accuracy of u n j. In all numerical examples, we 
observed that it is increasing in j, unless it reaches the peak level, which can be, in some cases, 
more or less the half of the arithmetic precision in a given computer algebra system. We used 
the system Maple 14 with 128-digit (decimal) arithmetic precision. 

Example 5.1. Consider the following CF expansion 



i> (-+3+,) + ^ (3*3=*) _ ^ _ ^ (s±l=H) 



x+K 



(2n - if 



x 



(2n) 



x 



(5.1) 

(3?cc > 0, cf. p. 33, Eq. (12)]), where ip(x) denotes the digamma function, tp(x) = V (x)/T(x). 
Let us put x := 1 and v := 1/2. We have the following CF of subclass 

{2n-l) 2 -\\. (2n) 2 



V := 1.327052799890558739735... = 1 + K 



n=l 



1 



1 



(5.2) 
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It is very slowly convergent. For example, the 100th classical approximant Sioo(0) = 1.319558 . . . 
yields <5ioo = 2.25. One can check that <5iooo = 3.24 and <5ioooo = 4.24, so it seems that 
accuracy of classical approximants is increasing logarithmically. Nevertheless, one can derive 
very good appr oxim ation of its value by using modified approximants S n (uj n ) with uj n « f^ n \ 
The Algorithm 4.1 provides the approximations u n j of tails of increasing accuracy m j, 

in the sense that 



UnO 
Unl 

U n 2 



T-2U, 



225 



r-2 11 + r„ + — -n- 1 +0(n 2 ) , 



T_ 2 n + T + T 2 n 1 + 



7125 
32768 



n 



-2 



868725 



1048576 



n 



-3 



+ 0(n" 4 ) 



u n3 = T_ 2 n + T + T 2 n + T^n + 



497625 
4194304 



n 



+ 0(n" 4 ) 



u 



+ 0[n-i), i>4, 

where r, are the coefficients of the asymptotic expansion of the tails u n = /C 2 "" 1 ) (see (|2j|). 



Let us observe that in each step of iteration in Algorithm 4.1, we derive only one (nontrivial) 
coefficient of this expansion; let us remark that Tn—\ vanishes for all i G Nq. This is the worst 



case to accelerate (cf. (3.11)). However, the approximation u^io yields about ten decimal digits 
of value V of fraction (5.2): 



Si (m m) 1.327052799 780862, and so 5 



1,10 



10.08. 



Table [3] shows the accuracies of all approximations u n j needed to be computed in order to 
obtain the value u% io- One can observe that accuracy of u± (see the second row of the table) 
is increasing linearly in j. Moreover, it is worth to analyse the Figure [TJ since it shows the 



Table 3: Values 5. 



derived by Algorithm 4.1 applied to CF (5.2). 

9 10 



n\j 





1 


2 


3 


4 


5 


6 


7 


8 


1 


1.24 


2.40 


3.24 


4.04 


4.82 


5.62 


6.44 


7.29 


8.17 


2 


1.79 


3.03 


3.90 


4.72 


5.53 


6.36 


7.22 


8.10 


9.03 


3 


2.13 


3.46 


4.38 


5.25 


6.11 


6.98 


7.87 


8.80 


9.78 


4 


2.38 


3.78 


4.76 


5.68 


6.58 


7.49 


8.43 


9.40 




5 


2.57 


4.04 


5.08 


6.04 


6.98 


7.94 


8.91 






6 


2.73 


4.25 


5.34 


6.34 


7.33 


8.32 








7 

8 


2.86 
2.98 


4.44 
4.60 


5.57 
5.77 


6.61 
6.85 


7.64 










9 


3.08 


4.74 


5.95 














10 


3.17 


4.87 
















11 


3.25 



















9.10 10.08 
10.01 



values 5 n j — accuracies of u n j — for large values of n. One can observe that each step of 
iteration in Algorithm 4.1 is of great importance, since it is much more productive than using 
approximations u n j for large values of n, provided that j is fixed. 



On the other hand, one can verify that even part of CF (5.1 ) is equivalent to the CF of the 
form K(a n /1), where 



a n + 



1 



1 



(x 2 



n 2 + 



1 



(3x 2 -3-^ 2 )n~ 3 + 0(n- 4 ), 



(5.3) 



4 16' ' 32 

and so it is of parabolic type. Thus, we cannot use 'improvement machine' of Lorentzen and 
Waadeland (see |j, |J p. 227]). Let us remark that 3f?x > 0, so CF is convergent by some 



14 



16 




Figure 1: Values 5 nj for CF (5.2) 



generalisation of Worpitzky's theorem (see (9j pp. 45-50], pEJ Thm. 3.30, p. 136]). Moreover, we 
get a n + \ = 0(n~ 3 ), since we took x = 1. Using the theorem of Thron and Waadeland [8], we 
show that modified approximants S n (—l/2) are convergent to the value V faster than classical 
approximants 5^(0), and so the 'fixed point' method works. However, it is not very efficient, 
since the accuracy of S n (— 1/2) behaves similarly to d n o for consecutive values of n. Namely, for 
n = 10,100,1000,10000, we get the following values of acc(5 ri (-l/2)): 4.52,6.36,8.35,10.34, 
respectively. In order to obtain 14 decimal digits of V we have to compute S67000o( — 1/2). This 
should be compared with <5i.i3 = 14.0 obtained after only 13 iterations in Algorithm 



4.1 



Example 5.2. Let us now set x := 1/2, v := 
0.883414269615 and even part being equivalent to CF of the form K(a n /1), where 



1/2 in (5.1). We obtain CF with value V 



a n + 



3 

64 



n- 2 + C(n" 3 ) 



(5.4) 



It is of parabolic type, too, hence we cannot use 'improvement machine' of Lorentzen and 
Waadeland. One can verify that conditions given by Thron and Waadeland in [H] are not 
satisfied, and we cannot use their results to accelerate the convergence, either. Indeed, 'fixed 
point' method fails — it gives only about 0.5 accuracy more than classical approximants S n (0). 
For example, we have S 10 = 0.78, acc(5i (-l/2)) = 1.18; 5ioo = 1-21, acc(5ioo(-l/2)) = 1.67; 
5iooo = 1-69, acc(5iooo(-l/2)) = 2.16; <5i ooo = 2.19, acc(5i ooo(-l/ 2 )) = 2.66. On the other 



hand, one can effectively approximate the value V of the considered CF using Algorithm 4.1 
The accuracies 5\j for j = 0, 1, . . . , 10 are: 



1.02, 2.15, 3.00, 3.82, 4.65, 5.51, 6.41, 7.33, 8.29, 9.29, 10.32. 



One can verify that CF in (5.1) is equivalent to a the fraction considered by Paszkowski [6]. 
Namely, we have 



x+ K 

n=l 



(2n 



+ 



(2n) 



~ x + - M 
x 



-1,-^,0,0 



(5.5) 
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where class Ai is defined in [61 Eq. (5.8)]. Using Paszkowski's analytical formulae gives a faster 
convergent CF. For example, its classical approximant /20 yields acc(/2o) = 21.6, while our 
method (applied to (5.1)) — #1^0 = 20.0. 

Example 5.3. Now, let us consider the following CF of subclass T>^ : 



f{z,a) := 



e v v a dv 



00 

z+K 

n=l 



n + a 



1 



n 



z > 0, a G R (5.6) 



-1/2 



(cf. [7J p. 103, Eq. (14)]). The asymptotic expansion of its tails contains powers of n 
(see Theorem 2.1), and thus we cannot use Paszkowski's method [6] in order to accelerate its 
convergence. Using its even part, we derive the following one- variant CF of subclass C2V 



f(z,a) 



a z 



z + l 



00 

!+Kr 

n=2 



-(n- l)(n + a- 1) 
2n + z + a — 1 



(5.7) 



Unfortunately, we cannot use the method of [5], since the above fraction does not belong to 
the subclass X\ C C21 (see the notation introduced in Eq. (2.3)]). Further, the equivalent 
transformation into limit periodic CF of the form K(a n /1) gives 



1 z 

"7 + 7 n 
4 4 



X + 0(n- 2 ). 



(5.8) 



Let us put z := 1/16 and a := 4. Then CF in (5.6) is convergent to 

V w 3.09147726049419952742569567195. 

We have £10 = 0.61, £50 = 3.02, 5ioo = 4.45 which shows quite slow convergence of the fraction. 
Applying 'fixed point' method to CF K(a n /1) equivalent to (5.7) seems to be worthless, since: 

6 W = 1.45, acc(Si (-f )) = 1.45, 
5 50 = 4.45, acc(5 50 (-^)) =4.53, 
<5ioo = 6.17, acc(S , iQ (-|)) = 6.23. 



Indeed, equation (5.8) does not follow the assumptions of Thron and Waadeland [5J. Let 
us remark that the computations of Algorithm 4.1 is organized so that array u n j is upper 
triangular. Thus, if we want to compute the approximation u± )r , then we need to start with 
r + 1 initial approximations u n o (n = 1, 2, . . . , r + 1). Applying our method to (5.6) yields uijg 
with accuracy 5\jq = 26.23. Let us note that S n (0) for n = 0, 1, . . . , 160 yield at most 6 exact 
significant decimal digits of V. 



Example 5.4. Now, let us consider CF of subclass 2?2l : 



log(l + x) 



1+K 

71=1 



2 

n x 
2n 



2 

n x 



2n+ 1 



x (— OO, —1] 



(5.9) 



(cf. [TJ p. 152, Eq. (7)]). We use its even part and obtain the CF of subclass C21, of the form 
K(c n /d n ) with 



22 ^2 ^2 ^2- 

-x n + 2x n x H — x n 

2 4 



O n 



d n = (2a + 4)n - 2 - x + (^y + 2^ n" 1 + 0(n~ 2 ) . 



(5.10) 
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One can easily check that it is of loxodromic type and belongs to subclass X\ (see notation in [5j 
Eq. (2.3)]). Thus, one can use either 'improvement machine' or the method described in [5]. Of 
course, one can also use 'fixed point' method, classical approximants of CF, as well. We try to 



compare their results with the results obtained by Algorithm 4.1 applied to CF (5.9). 
Let us put complex value of parameter x ■ 



-1.5 + O.Oli. Notice that it is very close to 



the ray (— oo, —1], which is the divergence region of (5.9). Hence, CF (5.9) is convergent very 
slowly. Indeed, we have 



'10 



-0.3, 5: 



50 



-0.5, 



noo 



0.1, 5. 



200 



0.6, 5, 



300 



0.9. 



We applied 'improvement machine' and the method described in |5] to CF K.(c n /d n ) (cf. ( 5.10| )), 
and the Algorithm 4.1 to CF (5.9). All of the methods construct upper triangular array of tails' 



approximations. We provide all of them with only fifteen initial approximations. Table [4] 
presents the obtained accuracies. One can observe that 'improvement machine' yields only 4 



decimal digits, while method of [5] and Algorithm 4.1 — 10 and 16 exact significant decimal 
digits of V, respectively. 



Appendix A Computing coefficients of asymptotic expansion of 
tails 



Here we show some formulae for coefficients r,- of the expansion (2.7) of tails u n for all considered 



subclasses of continued fractions (see Table [TJ . All the proofs will basically follow the idea given 
by Wynn in [10] . More precisely, we put the expansions (2.3) and (2.7) into the recurrence 



relation (2.5) in order to obtain the coefficients tj. After some rearrangement of terms, this 
leads to the following equation: 



oo Lf J+3 Lf J+2-i 

°= E E E liqjr m ~2i-2jn 

m=—8 i=—l j=— 1 

oo LfJ+2 



oo Lf J+ 4 m+4-2i 

m/2 + E E E^i 

m=— 10 i=— 1 j=— 4 



oo m+1 



+ Y E (Pi T m-2i - Pi^m-2i)n ™ /2 - } y p[ q m -j Tl "\ 

m=— 8 i=— 2 m=—3 i=—2 

where ojj are coefficients of the expansion of u n +i = ^°^_ 4 ooj n^^ 2 . One can check that 



(A.l) 



W_ 4 = T_ 4 , W_3 = T_ 3 , W_2 = 2t_ 4 + T_ 2 , ^0 = 7"-4 + T-2 + Tq, 





^ - 1 o + - Of) + - G"-D + • — ( _(2i : 1)/2 1 • i • 



= _\) + ti G -2) + ■ + n > ( ) ■ 1 e n - 

where (^) means generalised Newton's symbol 

r lvn (-z)(-« + l)---(-z + m-l) 



z£C,m GNU {0}. 
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Let c m denote the coefficients of powers of n m / 2 in (A.l) (m 
them are as follows: 



-10, —9, . . .); the first few of 



C-10 = Q-l T_ 4 , 

c_ 9 = 2(j_ x r_ 4 T_ 3 , 

c-8 = q'-i q-i t-4 + q'-i f_ 4 (2 r_ 4 + r_ 2 ) + 

+ q'-l T-3 + q'-l T -2 T_4 + % T 2 4 + T_ 4 (P-2 ~ P_ 2 ), 

C-7 = t_3 g_i + q_ 1 r_ 4 (t_i + 3r_ 3 /2) + q'_ 1 t_ 3 (2 r_ 4 + r_ 2 ) + 
+ r -2 q'-i r-3 + r_i g'_! r_ 4 + 2 r_ 4 g r_ 3 + r_ 3 p_ 2 - V-2 T -3- 



(A.2a) 
(A.2b) 

(A.2c) 
(A.2d) 



Equating them to zero gives the infinite system of equations satisfied by the unknown coefficients 



To: 



0, 



m 



-10,-9,..., 



(A.3) 



Taking into account the forms of numerators a n and denominators b n (see (2.3)), there easily 



follows that fi < 2 if CF0 T>2o- Hence, the first few equations of (A.3) are obviously satisfied. 



In Lemma A.l we analyse this in the case of all considered subclasses of CFs. In Lemma A. 2 
we show the form of quadratic equation satisfied by the beginning coefficient r^. 



Lemma A.l. The following formulae hold: 
1. For 



io- 



r_2 = 0, 



TO 



2p' ~ 2q' q +p- 1 - 2p 



2. In the case ofDf Q . 



T-l = T\ = T 3 



3. For CFeVn: 



TO 



P-l-P-1 
P-1<?0 

p'_l-p-l 



p-i+P^-po (p-i-P-i)<?i 
~r „' "T / 2 



T-i = ri = r 3 



TO 



'pU 

- qo 



«/r_ 2 = 0, 
i/r- 2 /0; 



p-i 



i/r_ 2 = 0, 
i/r_ 2 /0; 



4. For CFgVT?: 



5. For CF £ V. 



20- 



T-3 = T-l 



T-i = t_ 3 = 0; 



if t_ 4 = 0, 

T-2 = < 2p_ 2 +p'_ 1 -P-l . 3 i(p_2-p'_ 2 ) •/ / n 
H T2 if T-4 f= 0; 



1o 



6. For CF £ T>2i : 

T-3 = T-l = T\ = . . . = 0, 

p'-2Qo + p'-iq-i - {q'-iqo + q' q-i - + p-i - P-2) T -2 - (q'-i + q' ) T - 



TO 



2 4-\T-2 + P-2 - p'- 2 + q'-iq- 
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Proof. Each part of the lemma has to be considered separately. For the subclass T>io, the 
equations (A. 3) for m = — 10, —9, . . . , — 5 are trivially satisfied. For m = —4, we obtain that 



= t_2 {% r-2 ~ P-i + P-i) , 



so there must be r_2 = in the case of T> 1Q . Then the equation of (A. 3) with m = - 
evidently satisfied, and those corresponding to m = —2, —1 can be simplified to the form: 

<?o r -i - Q0P-1 = 0, 
r -i ( 2 Po - 4r o<?o - 2 QoQo ~ 2 Po+P-i) = 0. 



(A.4) 
—3 is 

(A.5) 



The latter equation yields the result of the part 1 . 

Using similar analysis, one can prove the part 4. Namely, in the case of class T>2o, the 
equations (A.3) with m = — 10 and m = —9 are obviously satisfied. However, equations 
corresponding to m = —8 and m = —7 can be written in the form 



t_ 4 (q' t_ 4 + P-2 ~ p'- 2 ) = 0, 
T-3 (2^o r_ 4 + P-2 - p'- 2 ) = 0. 



(A.6) 
(A.7) 



Hence, r_4 = in the case of T> 20 . Then equation ( |A.7[ ) is trivially satisfied, and the next 
equation of (A.3), related to m = —6, can be written as 



% T -3 



0. 



This proves the result of the part 4. 

The proofs of the remaining parts have a lot in common, so that we show only the proof of 
the part 2. Equation of (A.3), related to m = —3, can be simplified to the form 

r_i (2 %t_ 2 + P-i ~ P-i) = 0. 



Comparing this with (A.4) implies that r_i = 0, regardless of whether r_2 = or not. Formula 
for ro follows from the equation of the system (A.3) with m = —2. Next, using induction on 
j > 0, one verifies that the equation of (A.3), related to m = 2j — 3, can be simplified to the 
form 

T2j-\ (2<?o T -2 +P-i - P-i) = 0. 
Hence, all the coefficients T\, t%, . . . vanish. □ 



Lemma A. 2. The beginning coefficient r_o of expansion (2.7) satisfies the quadratic equation 



aT 2 + /3r + 7 = 0, 



(A.8) 



where a, /0, 7 are given in Table^ 

Proof. The proof uses Lemma A.l and the arguments in its proof. In the cases T>f and T>Yq the 
coefficients a,f3,j result from (A.4) and ( |A.5 ), respectively. Considering the subclass T>n, the 
equation (A.3) for m = —10,— 9, —8, —7 (cf. (A.2a)-(A.2d)) holds trivially, while for m = —6 
can be written in the form 

T-2 9-1 ( r -2 +q-i) = 0. 



In the case of subclass T> 2 q, using Lemma A.l we have that the equation of (A.3) with m 
is trivially satisfied. The next equation, related to m = —4, is of the form 



Qo T -2 + (p-i - P-i ~ P-2) r-2 - P-2 qo = 0. 
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In the case of £>20> the coefficients a, /3, 7 easily result from ( A.6). Finally, in the case of T>2i, the 
equations of the system (A. 3) related to m = —10, —9, —8, —7 (cf. ( A.2a)-( A. 2d)) are obviously 
satisfied, and the one corresponding to m = — 6 can be written in the form 



q'-i T-2 + (P-2 ~ P-2 + q'-i Q-l) T-2 ~ P-2 Q-l = 0. 



(A.9) 

□ 



The Lemma A. 3 extends the results given in the part 4 of Lemma A.l for the subclass T> 2 q. 
Notice that Lemma A. 3 follows by application of Theorem 2.1 which is proved using Lemma 

EH 

Lemma A. 3. In the case o/P^, we have 



T-l = Tl = T 3 



0. 



TO 



P-i go + P-2 gi ~ (gi + go) rl 2 + (p ; „ x + go - po - gp go) r_ 2 
2-r_ 2 g -P-i 



Proof. We continue the strategy of the previous proof. The equations of the system (A. 3) 
corresponding to m = —3 and m = —2 simplify to 



r_i ( 2t_2<7 + P-i ~ P-i ~ 2 p - 2 



0. 



and 

t (2 T_ 2 <7o + P-i ~ P-i) = P-iQo + P-2Q1 ~ tI^'q - (q[ + %) r 2 2 + (p'_ t + p - p - %q ) r_ 2 , 

respectively. The latter equation yields the formula for tq. One can check that the coefficient 
of r m in the equation of the system (A.3) related to m > —2 is given by 



2t_ 2 <7o + P-l ~ P-i + -jP-i- 



(A.10) 



Moreover, using the induction on j G No, one may verify that the equation of (A.3), corre- 
sponding to m = 2j — 1, can be simplified to the form 



, . . 2j - 1 

T2j-i 2r_ 2 g +P-1 -P-x H n — P~ 2 



0. 



To complete the proof, we show that expression (A. 10) cannot vanish for any integer m > —2. 
Namely, suppose that 

2t_ 2 q + P-i ~ P-i + -jP-2 = 



for some m £ N U {—1, 0}. Then we have 



Let us denote 



1 m _ 2t_2<? +P-1 - P-i 

2 " "¥ = V~2 

V / /3 2 -4a 7 



(A.ll) 



P-2 



where a, j3, 7 are given in the Table [2] (cf. the class ^20)- From Theorem 2.1, we derive that 
2r -2 % + P-i - p'-i = P-2 + sgn (5Rx) — 4aj 



(cf. (2.10)), and hence the inequality (A.ll) can be written in the form 



! TXX 

- > -— = l + sgn(5ftx) x. 



We conclude that x is a real number satisfying \ > 1 + |x| > 1, which is a contradiction. □ 
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Table 4: Computation of (5.9) with x = — 1.5 + O.Oli by the methods: 'improvement machine' 



[H p. 227], method of [5], and Algorithm 4.1 



(a) 'Improvement machine': values of acc(S n («;£ )) 



n\j 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


1 


0.6 


0.8 


1.1 


1.4 


1.6 


1.9 


2.1 


2.4 


2.6 


2.8 


3.1 


3.3 


3.6 


2 


1.4 


1.7 


2.1 


2.4 


2.7 


3.0 


3.3 


3.6 


3.9 


4.2 


4.5 


4.7 


5.0 


3 


1.8 


2.3 


2.7 


3.1 


3.5 


3.8 


4.2 


4.5 


4.8 


5.1 


5.4 


5.8 


6.1 


4 


2.0 


2.6 


3.1 


3.6 


4.0 


4.4 


4.8 


5.2 


5.6 


5.9 


6.2 


6.6 




5 


2.2 


2.9 


3.5 


4.0 


4.5 


4.9 


5.4 


5.8 


6.2 


6.6 


6.9 






6 


2.4 


3.1 


3.8 


4.3 


4.9 


5.4 


5.8 


6.3 


6.7 


7.1 








7 


2.5 


3.3 


4.0 


4.6 


5.2 


5.7 


6.2 


6.7 


7.2 










8 


2.7 


3.5 


4.2 


4.9 


5.5 


6.0 


6.6 


7.1 












9 


2.8 


3.6 


4.4 


5.1 


5.7 


6.3 


6.9 














10 


2.9 


3.8 


4.6 


5.3 


6.0 


6.6 
















11 


3.0 


3.9 


4.7 


5.5 


6.2 


















12 


3.0 


4.0 


4.9 


5.7 




















13 


3.1 


4.1 


5.0 






















14 


3.2 


4.2 
























15 


3.2 



























13 



14 



3.8 
5.3 



4.0 



n\j 





1 


2 


3 


(b) Method of 0: 
4 5 6 


values of acc (S„(t 
7 8 9 


*)). 
10 


11 


12 


13 


1 


0.6 


1.8 


2.5 


3.3 


4.0 


4.7 


5.4 


6.1 


6.8 


7.5 


8.1 


8.8 


9.4 


10.1 


2 


1.4 


2.4 


3.1 


3.9 


4.6 


5.3 


6.0 


6.7 


7.4 


8.1 


8.7 


9.4 


10.0 


10.7 


3 


1.8 


2.8 


3.6 


4.4 


5.1 


5.9 


6.6 


7.3 


7.9 


8.6 


9.3 


9.9 


10.6 




4 


2.0 


3.1 


4.0 


4.8 


5.6 


6.3 


7.0 


7.7 


8.4 


9.1 


9.8 


10.5 






5 


2.2 


3.3 


4.3 


5.1 


5.9 


6.7 


7.4 


8.1 


8.9 


9.6 


10.2 








6 


2.4 


3.5 


4.5 


5.4 


6.2 


7.0 


7.8 


8.5 


9.2 


10.0 










7 


2.5 


3.7 


4.8 


5.7 


6.5 


7.3 


8.1 


8.9 


9.6 












8 


2.7 


3.9 


5.0 


5.9 


6.8 


7.6 


8.4 


9.2 














9 


2.8 


4.1 


5.1 


6.1 


7.0 


7.9 


8.7 
















10 


2.9 


4.2 


5.3 


6.3 


7.2 


8.1 


















11 


3.0 


4.3 


5.5 


6.5 


7.4 




















12 


3.0 


4.4 


5.6 


6.7 






















13 


3.1 


4.5 


5.7 
























14 


3.2 


4.6 


























15 


3.2 





























14 



10.7 



n\j 





1 


2 


3 


4 


(c) Alt 
5 


;orithm 4.1 values of 5 n j- 
6 7 8 9 


10 


11 


12 


13 


1 


1.1 


2.4 


3.4 


5.0 


5.8 


6.9 


7.9 


9.0 


10.0 


11.1 


12.1 


13.1 


14.2 


15.2 


2 


1.8 


3.1 


4.4 


5.8 


6.7 


7.9 


8.9 


10.0 


11.0 


12.1 


13.1 


14.1 


15.2 


16.2 


3 


2.2 


3.6 


5.2 


6.4 


7.6 


8.7 


9.8 


10.9 


11.9 


13.0 


14.0 


15.1 


16.1 




4 


2.5 


4.0 


5.8 


6.9 


8.3 


9.4 


10.6 


11.6 


12.8 


13.8 


14.9 


16.0 






5 


2.7 


4.3 


6.3 


7.4 


8.9 


10.0 


11.3 


12.4 


13.5 


14.6 


15.7 








6 


2.9 


4.6 


6.7 


7.9 


9.4 


10.6 


11.9 


13.1 


14.2 


15.4 










7 


3.0 


4.8 


7.1 


8.3 


9.9 


11.1 


12.4 


13.7 


14.9 












8 


3.2 


5.0 


7.4 


8.7 


10.3 


11.6 


13.0 


14.2 














9 


3.3 


5.2 


7.7 


9.0 


10.7 


12.1 


13.4 
















10 


3.4 


5.4 


7.9 


9.3 


11.0 


12.5 


















11 


3.5 


5.6 


8.1 


9.6 


11.4 




















12 


3.6 


5.7 


8.4 


9.9 






















13 


3.7 


5.9 


8.5 
























14 


3.7 


6.0 


























15 


3.8 





























14 



16.2 
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